Asymmetric, helical and mirror-symmetric travelling waves in pipe flow 
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New families of three-dimensional nonlinear travelling waves are discovered in pipe flow. In con- 
trast to known waves (Faisst & Eckhardt Phys. Rev. Lett. 91, 224502 (2003), Wedin & Kerswell, 
J. Fluid Mech. 508, 333 (2004)), they possess no rotational symmetry and exist at much lower 
Reynolds numbers. Particularly striking is an 'asymmetric mode' which has one slow streak sand- 
wiched between two fast streaks located preferentially to one side of the pipe. This family originates 
in a pitchfork bifurcation from a mirror-symmetric travelling wave which can be traced down to 
a Reynolds number of 773. Helical and non-helical rotating waves are also found emphasizing the 
richness of phase space even at these very low Reynolds numbers. The delay in Reynolds number 
from when the laminar state ceases to be a global attractor to turbulent transition is then even 
larger than previously thought. 

PACS numbers: 47.20.Ft,47.27.Cn,47.27.nf 



Wall-bounded shear flows are of tremendous practi- 
cal importance yet their transition to turbulence is still 
poorly understood. The oldest and most famous ex- 
ample is the stability of flow along a straight pipe of 
circular cross-section first studied over 120 years ago 
[H. A steady, unidirectional, laminar solution - Hagen- 
Poiseuille flow [H IH - always exists but is only realised 
experimentally for lower flow rates (measured by the 
Reynolds number Re = UD/v^ where U is the mean ax- 
ial flow speed, D is the pipe diameter and v is the fluid's 
kinematic viscosity). At higher i?e, the fluid selects 
a state which is immediately spatially and temporally 
complex rather than adopting a sequence of intermedi- 
ate states of gradually decreasing symmetry. The ex- 
act transition Reynolds number Ret depends sensitively 
on the shape and amplitude of the disturbance present 
and therefore varies across experiments with quoted val- 
ues typically ranging from 2300 down to a more recent 
estimate of 1750 (0, i, @, 0, i, i, E, [HI, [H). A 
new direction in rationalising this abrupt transition re- 
volves around identifying alternative solutions (beyond 
the laminar state) to the governing Navier-Stokes equa- 
tions. These have only recently be found in the form of 
travelling waves (TWs) [13) [13] and all seem to be saddle 
points in phase space. They appear though saddle node 
bifurcations with the lowest found at Re = 1251. The 
delay before transition occurs {Ret > 1750) is attributed 
to the need for phase space to become sufficiently compli- 
cated (through the entanglement of stable and unstable 
manifolds of an increasing number of saddle points) to 
support turbulent trajectories. 

In this Letter, we present four new families of 'asym- 
metric', 'mirror-symmetric', helical and non-helical ro- 
tating TWs which have different structure to known so- 
lutions and exist at lower Reynolds numbers. The asym- 
metric family, which have one slow streak sandwiched 
between two fast streaks located preferentially to one 
side of the pipe, are particularly significant as they are 



the missing family of rotationally-asymmetric waves not 
found in fis', and have the structure preferred by the 
linear transient growth mechanism 15]. They bifurcate 
from a new mirror-symmetric family which can be traced 
down to a saddle node bifurcation at Re — 773. This 
figure substantially lowers the current best estimate of 
1251 for Reg - the Reynolds number at which the lam- 
inar state stops being a global attractor. The relative 
sizes of this new Reg and Ret for pipe flow are then 
more in line with plane Couette flow {Ret = 323 [l^ and 
Reo = 127.7 [13, [ll: Re based on half the velocity dif- 



ference and half the channel width) than plane Poiseuille 



flow {Ret « 1300 [ii,[23,|2l|,|2l|23 and Reg = 860 [Tl]: 



Re based on the mean flow rate and the channel width) . 
Beyond suggesting that Reg in plane Poiseuille flow can 
be significantly lowered, these latest discoveries highlight 
the substantial delay in Re between new solutions ap- 
pearing in phase space and the emergence of sustained 
turbulent trajectories. 

The new solutions were captured by inserting a fully 
3-dimensional spectral representation (Chebychev in s, 
Fourier in (p and z where (s, z) are the usual cylindri- 
cal coordinates aligned with the pipe) of the velocity and 
pressure field into the governing Navier-Stokes equations 
as viewed from an appropriately rotating and translating 
reference frame in which the TW is steady [l3]- The re- 
sultant nonlinear algebraic system was solved using the 
Newton- Raphson algorithm [3(|. To start the procedure 
off, an artificial body force was added to the Navier- 
Stokes equations (see [l^) designed to give streamwise- 
independent vortices and streaks of a finite amplitude. 
The size of the forcing was then adjusted to find a bi- 
furcation point at which the translational fiow symmetry 
along the pipe is broken. New finite-amplitude solutions 
to pipe fiow were found if this solution branch could be 
continued back to the zero- forcing limit. 



The TWs previously isolated [ij, |lJ] were induced us- 
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FIG. 3: Phase velocity C in units of [/ as a function of a 
for the mirror-symmetric modes (sohd hnes) and asymmet- 
ric modes (dashed) at 4 values of Re near the saddle node 
bifurcation at Re = 773. 



FIG. 1: Velocity fields for the asymmetric mode at Re = 2900 
(top) and the mirror-symmetric mode at Re = 1344 (bot- 
tom) (both at a = 0.75). An instantaneous state is shown 
on the left and a streamwise-averaged state on the right. 
The coloring indicates the downstream velocity relative to 
the parabolic laminar profile: red(dark) through white(light) 
represents slow through fast (with zero corresponding to the 
shading outside the pipe). In-plane velocity components are 
shown by vectors. The maximum and minimum streamwise 
velocities (with the laminar flow subtracted) and maximum 
in-plane speed for the asymmetric mode are 0.33, —0.42 and 
0.03 respectively while for the mirror-symmetric mode they 
are 0.31, -0.43 and 0.08 (all in units of U). 




FIG. 2: Instantaneous axial vorticity along one wavelength 
of the pipe for the asymmetric mode (top) and the mirror- 
symmetric mode (bottom) at the same Re and a as Fig 
[T] Contours are at -f/-60% of the maximum absolute value 
(green- light / red-dark) . 



ing a forcing that was rotationally symmetric under 

R„i : {u,v,w,p){s,4>,z) {u,v,w,p){s,4> + 2n/m,z) 

for some m = 2,3,4,5 or 6. As well as this rotational 
symmetry, all the TWs also possess the shift-&;-reflect 



symmetry 
S : {u,v,w,p){s,(l),z) 



{u, -V, w,p){s, —<p, z + T:/a) 



where a is the base axial wavenumber (so the periodic 
pipe is 27r/a long) and take the form u{s,4>,z,t) = 
tt(s, (j), z — Ct) where C is the a priori unknown axial 
phase speed of the wave. In contrast, new rotationally- 
asymmetric TWs were found by using a forcing function 
which created vortices with the radial velocity structure 
u{s,cl>) oc 5Re{e-i/^(l - s2)^^^i[l + cos(i^f )]e™'>} 
and hence distributed energy across a band of azimuthal 
wavenumbers. This choice led to a branch of asymmetric 
solutions whose component fast and slow streaks are pref- 
erentially located to one side of the pipe (see Figs [T] and 
[2]). These asymmetric TWs are S-symmetric and have 
one phase speed C along the pipe. They extend beyond 
Re = 5000 and originate at a pitchfork bifurcation (at 
Re = 1770 when a — 0.75) from a mirror-symmetric TW 
family (see Figs [T] and [2]) which satisfies the additional 
shift-&-rotate symmetry 

: (u, V, w,p){s, (j), z) {u, V, w,p){s, (j) + Tr,z + n/a) 

(coupled with the S-symmetry, this implies invariance 
under reflection in the line = ±7r/2). The mirror- 
symmetric solutions undergo a saddle node bifurcation 
at much lower Re: Re = 1167 at a — 0.75 going down 
to a minimum of Re — 773 at a — 1.44: see Fig. [31 
The friction factors associated with the upper branches 
of these new modes are much larger than the rotationally- 
symmetric modes 13, l3| (see Fig. |4]). 

Both new families pos sess the characteristic features 
of the TWs found in [isL [3| : essentially 2-dimensional 
fast streaks near the wall and aligned with the flow, slow 
streaks in the interior which are aligned on average with 
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FIG. 4: Friction factor A := 2D G/pU'^ against Re for the var- 
ious families of travelling waves where G is the mean pressure 
gradient along the pipe and p the density. The lower dashed 
line indicates the laminar value Ajam = 64/i?e and the upper 
dash-dot line indicates the log-law parametrization of experi- 
mental data 1/%/A = 2.0 log(i?e\/A). The labels are m values 
for the rotational symmetry Rm of the different TW families 
all drawn at the wavenumber which leads to the lowest saddle 
node bifurcation. The new TWs shown - mirror-symmetric 
modes (solid) and asymmetric modes (dashed) - correspond 
to m = 1 and a = 1.44. The bifurcation point is marked with 
a dot. The inset shows the phase velocity C (in units of U) 
versus Re for the two types of mode at a = 0.75 where the 
bifurcation was originally found and the optimum a = 1.44. 



the flow, and a smaller (typically by an order of mag- 
nitude) 3-dimensional wave field. By continuity, helical 
TWs should exist with these fast streaks inclined to the 
flow direction and indeed a surface of such solutions can 
be found connecting the upper and lower branches of the 
mirror-symmetric TWs (see Fig. [5]). These helical TWs 
take the form m(s, t) = u{s, (p — j3[z — Ct] —ut, z — Ct) 
with /3 measuring the helicity in the Galilean frame mov- 
ing at Cz and w being an azimuthal phase speed relative 
to the Galilean frame. Helicity destroys S-symmetry but 
a modified form of 0-symmetry (0^) is preserved where 
the rotation transformation is now — > 4) + {\ — §)7r. 
The helicity (3 and rotational speed lo never rise above 
O(10~^) for Re < 1500 confirming the flow preference 
for non-rotating, axially-aligned streaks. Interestingly, 
in the range Re = 1165 — 1330, the helicity (3 on this sur- 
face passes through zero twice in going between the two 
mirror-symmetric branches (see Fig. [5]). These points 
correspond to an isola in the (fixed a) C vs Re plane 
of rotating non-helical modes which are neither shift-&- 
reflect symmetric nor have any rotational symmetry. The 
helical and non-helical rotating waves look very similar to 
the mirror-symmetric modes except for a slight twist in 
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FIG. 5: A schematic picture of how all the new travelling 
wave branches fit together in {(3,Re,C) space (at a — 0.75). 
The main parabolic curve in the /3 = plane is the mirror- 
symmetric branch off which the asymmetric branch bifur- 
cates (uppermost line). Helical branches bulge out of the 
P = plane and connect upper and lower parts of the mirror- 
symmetric. Across a finite range of Re, these helical modes 
perforate the /3 = plane in between the mirror-symmetric 
branches creating an isola of non-helical rotating TWs (closed 
dash-dot loop). Helical waves also connect the asymmet- 
ric branch and the helical solutions which originate from the 
mirror-symmetric solutions. 




FIG. 6: Two velocity slices across a helical mode taken at the 
same instant of time but 25 D apart with a — 0.75, P — 0.019, 
uj = —0.0011 at Re — 1344. The velocity representation is as 
in Fig. □ 



the streak structure along the pipe (see Figs [5] and [7]). He- 
lical modes continued off the asymmetric modes have no 
symmetry at all and originate in a symmetry-breaking bi- 
furcation off the 0/3-symmetric helical solutions extended 
from the mirror-symmetric waves: see Fig. [5l 

The asymmetric, mirror-symmetric and helical TWs 
all represent saddle points in phase space with a very 
low-dimensional unstable manifolds (e.g. 2 for the asym- 
metric mode at (a. Re) = (0.75, 1820) and 4 for the 
mirror-symmetric mode at (a. Re) = (0.75, 1184)). Their 
presence indicates the richness of phase space even at 



FIG. 7: The four fast streaks of the hehcal mode shown in Fig. [6] plotted over one (5 wavelength ~ 170 D. 



Reynolds numbers approaching 773. The delay of transi- 
tion until Re > 1750 suggests that the establishment of a 
'turbulence-bearing' scaffold constituted of all their sta- 
ble and unstable manifolds is far from immediate. The 
clear implication is that while the emergence of alter- 
native solutions to the laminar state seems a necessary 
precursor for transition, it is not a good predictor of the 
actual Reynolds number at which this occurs in pipe flow 
(and other shear flow systems). Once the transitional 
regime has been reached however, there is now mounting 
experimental [2^ and numerical evidence [2^, in- 
dicating that some of the travelling waves at least (those 
with low to intermediate wall shear stresses 26|) appear 
as transient but recurrent coherent structures within the 
flow. Intriguingiy, numerical simulations 26] have also 
revealed that a number of travelling waves with low wall 
shear stress sit on a dividing surface in phase space (a 
separatrix if the turbulence is a sustained state) which 
separates initial conditions which directly relaminarise 
and those which lead to a turbulent episode. Recent com- 
putations 2^, 2^ using a shooting technique to converge 
onto this dividing surface appear to have already found 
that the asymmetric wave sits there too (compare Fig 
[T]to Fig. 8 of [iil and Fig. 1 of [28|). The fact that 
this wave bears some resemblence to m = 1 'optimal' 
disturbances which emerge from linear transient growth 
analyses also suggests an enticing opportunity to bridge 
the gap between linear and nonlinear approaches. 

In summary, we have presented a series of new travel- 
ling wave solutions to the pipe flow problem which have 
different structure to existing solutions and which exist 
at far lower Reynolds numbers. One type - the asymmet- 
ric modes - represents the missing m — 1 family from the 
waves found initially These waves also appear to 

rationalise some interesting results from recent numerical 
computations ^2^, i22j , thereby corroborating the picture 
which is emerging that lower branch TWs (and therefore 
also their stable manifolds) sit on the separatrix between 
laminar and turbulent states. 

We acknowledge encouraging discussions with Fabian 
Waleffe and the support of EPSRC. 
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